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We study an elasticity model for compressed protein monolayers or particle rafts at a liquid 
interface. Based on the microscopic view of hard-core particles with soft shells, a bead-spring model 
is formulated and analyzed in terms of continuum elasticity theory. The theory can be applied, for 
example, to hydrophobin-coated air-water interfaces or, more generally, to liquid interfaces coated 
with an adsorbed monolayer of interacting hard-core particles. We derive constitutive relations for 
such particle rafts and describe the buckling of compressed planar liquid interfaces as well as their 
apparent Poisson ratio. We also use the constitutive relations to obtain shape equations for pendant 
or buoyant capsules attached to a capillary, and to compute deflated shapes of such capsules. A 
comparison with capsules obeying the usual Hookean elasticity (without hard cores) reveals that 
the hard cores trigger capsule wrinkling. Furthermore, it is shown that a shape analysis of deflated 
capsules with hard-core/soft-shcll elasticity gives apparent elastic moduli which can be much higher 
than the original values if Hookean elasticity is assumed. 


I. INTRODUCTION 

Many soft matter systems exhibit elastic properties that 
go beyond the Hookean linear elasticity. There are a num¬ 
ber of prominent examples of elastic materials with unique 
properties, which require tailored elasticity models for an 
accurate theoretical description. An early example is the 
Mooney-Rivlin law for large deformations of incompress¬ 
ible, rubber-like materials |T. To describe the bending of 
lipid membranes, the Helfrich energy was introduced [2] , 
and Skalak et al. and Evans developed a strain-energy 
function for deformed red blood cell membranes [M|. 
Only the combination of both types of elastic energy 
functionals with a Helfrich energy describing the lipid 
membrane and a Skalak strain-energy function describing 
the spectrin skeleton is able to describe the experimentally 
observed sequences of red blood cell shapes successfully 
111 El- The correct sequence of shapes cannot be repro¬ 
duced using a simpler effective model, for example, of 
the Helfrich type by choosing effective bending moduli 
properly. 

In general, indications whether simple elastic models 
such as a linear elasticity are sufficient to describe the 
deformation of a certain material deformation can be 
found when comparing experimental results to theoretical 
predictions. If, for example, linear elasticity is assumed 
in the theoretical modeling but comparison with experi¬ 
mental data shows that the resulting elastic moduli are 
apparently changing throughout the deformation, this 
suggests that a more advanced and detailed elastic model 
should be used. 

In this Article, we develop an elasticity model tailored 
for monolayers of particles or molecules at a fluid interface 
under compression. Since the pioneering work of Pickering 
[8], it is known that adsorbed particles at a liquid interface 
act as surfactants and can stabilize droplets in emulsions 
depending on the wettability properties [9]. Adsorbed 
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particles also stabilize foams nni Colloidal particles 
adsorbing at a liquid droplet interface form colloidosomes 
El, where particles arrange to spherical colloidal crystals 
[12] . or “armoured droplets” if particles are jammed m- 
A liquid interface coated with adsorbed particles starts 
to exhibit elastic properties typical for a two-dimensional 
solid or elastic membrane, such as resistance to shear, 
buckling, wrinkling and crumpling mm- This has also 
been reported for protein-coated bubbles [15], droplets 
m or vesicles [T8j. 

Hydrophobins are small proteins from fungal origin 
m , which adsorb strongly to the interface of an aque¬ 
ous solution because of their amphipathic nature with 
a hydrophobic patch on their surface [2U] , Layers of hy- 
drophobin can be used very efficiently to stabilize bubbles 
and foams [2H423] . which makes them interesting as a 
model system for protein particle rafts at the air-water 
interface and for applications. 

Particles coating liquid droplets have applications in 
food industry [MJ. Because of their strong amphipathic 
nature and biocompatibility, hydrophobins coating air- 
water interfaces have various applicationsJUJEIiES] rang¬ 
ing from medical and technical coatings [25] to the produc¬ 
tion of protein glue and cosmetics, or in the food industry 
in the stabilization of emulsion bubbles EZI- 

In a recently published experiment 128) , the deformation 
of a hydrophobin-coated bubble rising from a capillary was 
investigated. As a result of deflating the bubble, wrinkles 
appeared on the bubble surface proving that the protein 
layer has elastic properties and a nonvanishing shear mod¬ 
ulus. Furthermore, the analysis of this experiment showed 
that the elastic response is very nonlinear, and a steep 
increase in the measured elastic modulus with increasing 
bubble deflation was observed. It was conjectured that 
this steep increase is related to the molecular structure of 
the hydrophobin proteins consisting of a hard core (a /3 
barrel) with a softer shell (loop and coil structures). The 
hard cores coming into contact could then be the cause 
for the steep increase in the elastic modulus. 

In order to verify this interpretation, we develop an 
elasticity model based on the microscopic view of a particle 
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of such compressed particle-decorated liquid interfaces 
therefore requires a model as presented here, which prop¬ 
erly connects Hookean elasticity of soft springs or soft 
particle interactions with hard-core elasticity, and goes 
beyond a simple effective description based on Hookean 
elasticity. 


FIG. 1. Bead-spring model with hard cores for interacting 
hard particles at the air/water interface or a liquid interface. 
The springs give the network a Hookean elasticity, and the 
hard cores (green disks, diameter L ) impose limits on the 
maximal admissible compression (see configuration on the 
right). The lattice constant is normalized to 1 and L < 1. 


raft sketched in Figure [l] Globular particles interact by 
soft springs (corresponding to an outer soft shell) and 
steric interactions (corresponding to hard cores). From 
this microscopic view, we derive continuum elastic laws, 
which can then be used to analyze the deflated shapes 
of capsules with this type of elasticity. The hypothesis 
underlying this approach is the following: The steep 
increase in the surface Young or stretch modulus that was 
obtained in ref [28] is a result of applying linear Hookean 
elasticity in a situation where the proposed elasticity 
model for the protein particle raft that includes hard 
cores should be more appropriate. Therefore, we should 
be able to describe the observed shapes of hydrophobin- 
coated bubbles more exactly and with much less variation 
of the elastic moduli along the deflation trajectory if the 
new hard-core/soft-shell elasticity model is used. 

The proposed elasticity model is quite generic: It is 
a generalization of a standard bead-spring model for 
Hookean elastic membranes [29J, which also takes into 
account hard cores of the beads. On the other hand, 
it generalizes a particle raft elasticity model presented 
in ref [IT: which considers exclusively hard-core interac¬ 
tions. Therefore, our hard-core/soft-shell elasticity model 
applies more widely than to monolayers of the protein 
hydrophobin. We expect such contact interaction to be 
relevant for many liquid interfaces decorated with inter¬ 
acting hard particles [13], ranging from colloidal rafts 
covering colloidosomes EJ or “armoured droplets” [13] 
to rafts of larger particles [14] . 

Our results will show that the elastic properties of the 
two-dimensional material are Hookean for small compres¬ 
sion, where only soft springs are loaded without hard-core 
contact but strongly modify as soon as compression brings 
the hard cores of particles into contact. The constitutive 
stress-strain relations of the soft spring, which determine 
the material properties for small compression, have to be 
properly connected to stress relations, which are governed 
by force balance for compressions that induce hard-core 
contact. The transition between both regimes can take 
place between different regions within the same material, 
as we will see in the analysis of shapes of pendant or 
buoyant capsules attached to a capillary. A detailed and 
quantitative understanding of shapes and deformations 


II. CONTINUUM DESCRIPTION OF A 
BEAD-SPRING MODEL WITH HARD CORES 


We assume that the beads are arranged in the (x, y)- 
plane in a hexagonal crystal with a lattice constant nor¬ 
malized to 1, see Figure [T] Such an arrangement is the 
closest packing of spheres and it behaves isotropically 
[3D], so that the continuum elasticity of the membrane 
can be characterized by two elastic moduli, for example 
the surface Poisson ratio v and surface Young modulus 
F/ 2 D ■ On the microscopic level, the elastic response of 
the membrane is governed by the spring constant k. By 
evaluation of the deformation energy within one unit cell, 
it can be shown that the continuum elastic moduli are 
determined by jDDJ f5Tj 



and 



(1) 


Without the hard-core interactions, the membrane can 
thus be described with the usual Hookean elasticity as 
specified by the strain-energy density [32] 


Ws(\x,\y) = l^^[(\ x -1) 2 

+ 2v(X X — l)(Ay — 1) + (Ay — l) 2 ] + X X Xy"f, 


where X x and X v denote the stretches in x- and y-direction, 
respectively. The strain energy density is measured per 
surface area of the undeformed interface. The form of the 
surface energy density in eq[ 2 ]is valid only for deformations 
where the x- and y-directions are the principal directions 
of the strain tensor, which we assume throughout this pa¬ 
per. This means that a line element (da:o,d?/o) is mapped 
to (dx,dy) = (A x da;o, Aydj/o) by the deformation. 

The term XxXy'y accounts for an isotropic surface ten¬ 
sion 7 acting in the liquid surface [28] . We assume that the 
corresponding liquid interface covers the entire deformed 
area X x X v Aq, where Aq is the undeformed reference in¬ 
terface area. In the presence of adsorbed particles, which 
are only partially wet, a fixed area Ah ar d of the liquid 
surface is replaced by a liquid-solid-interface with a sur¬ 
face tension 7 hard, and the last term in eq [ 2 ] becomes 
{X x Xy A h ar d j A q Yf M |, //loh'Uird ■ This shift by a 
constant independent of the stretching factors does nei¬ 
ther change the interface stresses nor the resulting shape 
equations. 

From the energy density eq[2j the stresses acting in the 
interface can be derived [53] , A superscript ^ indicates 
that this is the contribution of the springs, 

ri s) (A„ Ay) = [(A x - 1) + v{X y - 1)] + 7 , (3) 

JL V Ay 
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FIG. 2. Forbidden domains (light gray) in the (A^, Aj,)-plane limited by ellipses, (a) Different orientations ipo produce different 
forbidden domains. Straight lines are degenerate ellipses with infinite major axis, and dashed lines with alternating color 
indicate if two of the three ellipses are identical. The pictograms show the lattice orientation, (b) Choice of ipo that minimizes 
the area of the forbidden domain. For easy identification, the boundaries are labeled 1 and 2; and the point on both lines is 
termed B. There is an additional vertical and horizontal boundary at A i/L = l/%/3 where next nearest beads come into contact, 
(c) Sketches of the different possibilities of contact between the beads, here for tpo = 30°, that is, when the compression is 
predominantly in ^-direction. 


and with indices x and y interchanged for t^ s \ We note 
that compression with Ai < 1 gives rise to negative elastic 
contributions to the stresses, whereas the surface tension 7 
always provides a positive contribution, such that t* < 7. 

Now we have to evaluate how these results from the 
spring elasticity are modified by the steric interactions 
between the hard cores. The springs in the lattice have a 
rest length of 1 and can be oriented along three different 
directions i, which we characterize by the angle (fi to the 
i-axis. This angle can take the values 77 = </?o+*tt/ 3, with 
i £ {— 1 , 0 , 1 }, where ipo determines the overall orientation 
of the lattice in the {x, y)-plane. Upon deformation, the 
length of a spring along direction i changes from 1 to 


/ X x cos tpA 
\Xy sin ipi J 



The steric interactions enforce that the springs of the 
lattice can be compressed at maximum to a minimal 
length of L (with L < 1), which corresponds to the 
diameter of the hard cores, measured in units of the 
lattice constant. Thus, we have three conditions di> L 
to be satisfied, or equivalently 

A 2 cos 2 (v?o — 7t/3) + A 2 sin 2 (vj 0 — 7t/3) > L 2 

A 2 cos 2 (v? 0 ) + A 2 sin 2 (<^o) > L 2 (5) 

A 2 cos 2 (<£o + tt/3) + A 2 sin 2 (v? 0 + 7 r/3) > L 2 . 

In the (Aa;, A !/ )-plane, these conditions specify three el¬ 
lipses to be excluded from the admissible domain for the 
stretches; see Figure [2jr. All three ellipses contain the 
point (L, L ). 

The excluded or “forbidden” domain of stretches, 
shaded light gray in Figure E* and b, depends on the 
orientation <po of the lattice. We choose this parameter 
according to the following rule: If \ x < X y , then ip 0 = 30°; 
otherwise ipo = 0°. With this choice, the forbidden do¬ 
main becomes as small as possible; see Figure [2]:>. Choos¬ 
ing ipo in dependence of the stress state means that it may 
change during a deformation; for particle rafts this seems 
plausible because the particles are not rigidly cross-linked, 


but may rearrange to change their lattice orientation from 
a 0 ° to a 30° state (see pictograms in Figure [ 2 J;). If the 
stress state of the membrane becomes inhomogeneous, 
as it will happen, for example, for capsules attached to 
a capillary as considered below in more detail, we may 
encounter regions of 0° and 30° orientation within a single 
membrane. Then, problems might arise at the boundary 
separating two regions of different orientation, because 
the lattices cannot be joined properly. This could give 
rise to a line energy; we neglect these complications in 
the following. 

With our choice of ipo, the boundary line between the 
forbidden and admissible domain of the stretch plane can 
be described by (AfA^A^) = (X^(X y ), X y ) with 


Ai b) (A y ) 


yj4L*/3->%/3 

t/^SA 2 


if A y > L (boundary 1) 
if A y < L (boundary 2). 


( 6 ) 

Here we have introduced the terms boundary 1 and 2, 
which have to be distinguished in most of the following 
calculations. They are plotted in red and blue, respec¬ 
tively, in Figure [Ijr. The point (X x ,X y ) = {L,L) which is 
on both boundary lines is termed point B. 


If the external loads try to push the membrane into 
the forbidden domain, the hard-core interactions keep it 
on the boundary by providing additional contributions 
to the stresses t x and t v . These hard-core contributions 
are denoted by rjf) and The complete stresses then 
read 


t x = ri s) (X x , X y ) + r| c) and t v = r^ s) (A x , X v ) + r^ c) 

(7) 

with ( X x ,X y ) being a point on the boundary of the for¬ 
bidden domain as given by eqlbl Since and are 
transmitted through the “skeleton” of hard cores, they 
must satisfy certain conditions of force equilibrium that 
can be derived from the geometry of the lattice. In Ap¬ 
pendix |Aj it is shown that force balance prescribes the 
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FIG. 3. Admissible domains (light green) and hard-core 
interaction domains for stretches (a, c) and stresses (b, d) for 
L = 0.95 and i? 2 D = 57- The top row shows large views of 
the stress and stretch planes, and the bottom row close-ups 
of the relevant region. In (b) and (d), the straight thin lines 
indicate the stresses that can be reached from a point on the 
boundary by adding the hard-core contributions to the elastic 
stresses. The wrinkling region is shaded in gray in (d). The 
arrow in (c) indicates the trajectory of uniaxial compression 
of a planar layer considered below. 


On boundaries 1 and 2, the hard cores come into contact. 
In the stretch plane, this boundary cannot be trespassed 
because of the geometric constraints imposed by the un- 
deformable hard cores: Even if the external forces try to 
push the lattice beyond this line, the lattice will get stuck 
on the boundary of the forbidden domain. In the stress 
plane, however, the points beyond boundaries 1 and 2 can 
be accessed by including the hard-core contributions r| c ) 
and rin the stresses because hard cores can transmit 
force even though they are undeformable. 

A point (A^A^) on the boundary in the stretch 
space is mapped by Hooke’s law to a point ) 

on the boundary in the stress space. From this point 
on, stresses (■ t x ,t v ) = (r^,r^) + (t^ c \t^) can be 
reached, where the hard-core contributions r( c ) and r 
must be negative (because the skeleton can only support 
compressive stresses) and must obey the force balance 
constraint (eq [ 8 ]) . The submanifold of stresses which 
is accessible from the point on the boundary is, thus, a 
straight line with a slope X y /3X x starting from boundary 1 
or 3X y /X x starting from boundary 2 in Figure [3|i (see the 
thin red and blue lines). The analogous discussion of the 
boundaries where next nearest beads come into contact 
is spared because it is not relevant for the applications 
presented below. 

For point B, with X x = X y = L, the ratio of the hard¬ 
core stresses is not fixed to a certain value. Instead, it 
can range from 1/3 to 3. In the plot of the stress plane, 
Figure Si. this means that the whole white region is 
accessible from point B. 


ratio of the hard-core stresses to 

T y C ' > f A^/3 X x if X y > L (boundary 1 ) 
r i c) jsA^/A^ if Ay < L (boundary 2) 

1 1 J 
and - < —/y) — 3 if A x = A y = L (point B). 

3 Tc 

Note that (X x , X y ) must always be a point on the boundary 
of the forbidden domain when these equations are applied 
to calculate the hard-core contributions to the stresses. 

For strong compressive strains with A^, < l/y/3 or 
Ay < 1 /\/3, next nearest beads, which are not connected 
by springs, come into contact; see Figure [2jo, orange and 
violet lines. Such compressive strains are not relevant for 
the applications discussed in this Article. 

The above results for our custom elasticity model for 
particle rafts are summarized and illustrated in Figure [3] 
Figure^ and b shows a large view of the stretch plane 
and stress plane with different regions. Figure [3jc and d 
focuses on the regime relevant for compressed membranes: 
A i < 1 and 0 < n < 7 . The light green regions are 
the admissible regions, where the hard cores are not in 
contact. Here, the usual Hookean elasticity (eq[3]) is valid, 
and there is a bijective mapping between (A^, X y ) points 
in the stretch plane and ( t x ,t v ) points in the stress plane. 


III. COMPRESSION OF PLANAR FILMS 

To get a first impression of the elastic behaviour of 
a particle layer with hard cores at a fluid interface, we 
investigate a typical compression experiment in a Lang¬ 
muir trough [33 1 134 ] . In this setup, a monolayer is spread 
on a water surface and subsequently it is compressed 
in x-direction by moving the barriers of the trough. A 
Wilhelmy plate can then measure the stresses t x and t v . 

If we consider X y = 1 as fixed and compress the layer 
with a ratio A^ < 1 , we expect the monolayer to wrinkle 
at sufficiently high compression. The determination of the 
critical compression is a purely geometric problem. In the 
stretch space of Figure^, we follow a horizontal line until 
we cross boundary 1 as indicated by the arrow. Then 
the layer must wrinkle because the hard cores cannot be 
compressed. 

On the other hand, if we consider the compressive stress 
t x < 0 to be given, we can determine the critical stress 
from the stability equation mm 

E B d x w(x) — r x d x w{x) + pgw( x) = 0, (9) 

where E B is the bending stiffness of the layer, w( 2 ) is 
the displacement in the ^-direction which is assumed to 
be independent of y, and pg is the fluid density times 
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acceleration of gravity. The bending rigidity E B can 
arise from an additional bending stiffness of connecting 
springs. For hydrophobin layers this corresponds to a 
bending rigidity of the soft shell of the protein. Assuming 
sinusoidal wrinkles of the form w ~ sin kx, eq [9] gives a 
critical stress t x which still depends on the wavenumber 
k of our ansatz. Minimizing the magnitude of the stress 
with respect to k gives a critical wavelength and critical 
stress of EH 

A c = 2n ( E B /pg ) 1/4 and r x , c = -2 ^E B pg. (10) 


1.4 
1.2 
1.0 
a o.8 
^ 0.6 
0.4 
0.2 
0.0 

-1.4 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 

t x /E 2 o 



This is the same result as that for an elastic membrane 
without hard cores. The result for A c differs from what has 
been obtained in ref Il4t we do not find a dependence on 
the packing fraction of hard particles because the bending 
modulus corresponding to the soft springs appears in eq[T 0 | 
rather than an effective modulus depending on packing 
fraction as in ref HU 

Differences between membranes with and without hard 
cores can be found when the Poisson ratio of the layer 
is measured. For uniaxial compressive strain in the re¬ 
direction (A x < 1, A y = 1), the Poisson ratio of a linear 
elastic material can be measured as u = t v It x . Now let 
us apply this “measurement” to the more complex total 
stresses in eq[7j 


r(s) 


r( c ) — - 


T' + T » 7 

T i S ^ + T x C ' > ~ 7 


* I y 
IT. 


( 11 ) 


Here, we subtract the fluid surface tension 7 from the 
total stresses, which provides a background tension and 
is typically subtracted in measurements of the surface 
pressures Tl x and Ii y [33) [34]. For simplicity, we also 
neglect the 1/A X prefactor in relation (l3j) for t^ s \ which 
represents a geometrical nonlinearity. In The regime where 
the hard cores are in contact, the apparent Poisson ratio 
thus reads 


^(AW-^ + fW^AW ) 2 
i-1) + ^l c) 


( 12 ) 


with A), b ) = ^/4L 2 /3 — 1/3 and = t^/E 2 b- The 
hard cores have a nonvanishing contribution only if 
the applied compressive stress satisfies t x < (A^N — 1) 
so that the hard cores are in contact. For larger t x or weak 
compression, the usual Hookean law is valid, which leads 
to z/ app = v = 1/3. In the limit of strong compression 
where the hard-core contribution dominates the stresses, 
we have ^ app = l/3(A£ b l ) 2 = 1/(4L 2 — 1 ), which is always 
larger than 1/3 for L < 1 and can even become larger 
than 1. For L < l/-\/3, this relation changes because 
next nearest beads come into contact. This shows that 
contact of hard cores can give rise to a more effective 
redirection of compressive stress into the perpendicular 
direction, which gives rise to a stress-dependent Poisson 
ratio increasing with compression. The crossover between 
weak compression with the Hookean value v = 1/3 and 
the increasing apparent Poisson ratio r/ app for strong com¬ 
pression is shown in Figure[4]for different values of L. Our 


FIG. 4. Apparent Poisson ratio u app = H. y /H x in our elasticity 
model for different values of L as indicated. At the kinks, the 
hard cores come into contact and the hard-core contributions 
to the stresses begin to influence this “measured” Poisson 
ratio. 


results differ from i/ app = l/-\/3, which has been found in 
ref HU for pure hard-core particle rafts. In refGjU on the 
other hand, it has been argued that ^ app = 1/3 is correct 
also for pure hard-core particle rafts. Our result (eq[l2| 
for the more general elasticity model explains that, for 
rafts of interacting hard particles, where the interaction 
provides a “soft shell”, experimental measurements such 
as those in refs [331 and [33 should be interpreted using a 
stress-dependent Poisson ratio, which increases beyond 
1/3 and up to v app = 1/(4 L 2 — 1) if hard cores come into 
contact; see Figure [4] 

In ref 1371 for example, an experiment is reported where 
a monolayer of hydrophobin molecules is compressed in 
a Langmuir trough. Two Wilhelmy plates were used to 
measure the stresses in x- and y-directions, respectively, 
and from the data one can calculate an apparent Poisson 
ratio of zz app ss 0.6 — 0.7. This relatively large value can 
be well explained with our elasticity model if L < 0.8 
according to Figure [4] It is tempting to compare this 
hard-core length to available molecular information. In 
ref [3U it has been observed that visual buckling, which 
could be due to hard-core contact, happens at a molecular 
area of abuckie ~ 347 A per HFBII hydrophobin protein. 
Assuming closed-packed hard cores, this corresponds to a 
diameter of 20 A, if HFBII is assumed to be spherical. A 
value L = 0.8 then corresponds to a total diameter of 25 A 
of the protein. This value is consistent, for example, with 
the dimensions 24 x 27 x 30 A given for the whole HFBII 
molecule as obtained from X-ray crystallography [23] , 
Grazing-incidence X-ray diffraction and reflectivity results 
in ref[35]hint, however, at smaller diameters 20 —24 A and, 
thus, larger values of L. In refs (571 and f39l plateaus in the 
surface pressure isotherms have been observed, which can 
be interpreted as a liquid-gas coexistence of a gas phase 
of isolated hydrophobins and a dilute liquid where soft 
shells of hydrophobins just start to make contact. Taking 
the area per molecule a p i ateau at the high-density end of 
these plateaus we find values of L = (a p i ate au/“buckle ) 1 ^ 2 
in the range L = 0.85 — 0.90. Overall, a reliable estimate 
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capsule shell and are given by [58] 

r'(s 0 ) = A s cosip 
z'(sq) = A s sin^ 

$'{sq ) = -(p-Apgz-K v ,T ip ) ( 13 ) 

T s 

t'(s 0 ) = A s———cos ip. 


FIG. 5. Arc-length parametrizations in cylindrical coordinates 
of the undeformed (ro(so), «o(so)), deformed (r(s),«(s)) and 
the wrinkled midsurface. The wrinkled region of length L w 
is described by an axisymmetric pseudosurface around which 
the real midsurface oscillates. 


of L based on molecular information appears difficult. 


IV. SHAPE EQUATIONS IN THE 
PENDANT/RISING CAPSULE GEOMETRY 


We now apply the previously developed elasticity model 
to the shape equations for axisymmetric capsules attached 
to a capillary which have been developed in ref HU This 
is a challenging problem because it turns out that the 
transition between the elasticity regimes A, 1, 2, and B 
can take place within the same material between different 
axisymmetric regions along the capsule. Therefore, we 
have to generalize the shape equations from ref HU to 
include switching between different constitutive elastic 
laws corresponding to the Hookean regime A and the hard¬ 
core dominated regimes 1, 2, and B along the capsule 
contour. 

The shape equations are derived from non-linear shell 
theory and their solution describes the shape and stress 
distribution of a deformed pendant or rising capsule. Fig¬ 
ure [5] shows the parametrization of the undeformed and 
deformed capsule shape. At its upper rim, the capsule is 
attached to a capillary of diameter a. The reference shape 
(see Figure [5^,) is a solution of the Laplace-Young equa¬ 
tion with an isotropic interfacial tension 7 and a density 
difference A p of inner and outer fluid [58] ■ This shape 
can be deformed, for example by reducing the internal 
capsule volume. Each point (r’o(so), -o(so)) is mapped 
onto a point (r(so), z i s o)) in the deformed configuration, 
which induces a meridional and circumferential stretch, 
A s = ds/dso and Ap = r/ro, respectively. The arc length 
element ds of the deformed configuration is defined by 
ds 2 = (r'(s 0 ) 2 + /(so) 2 )dsg. 

The shape equations that determine the deformed con¬ 
figuration describe force and torque balance within the 


Here, ip denotes the slope angle as defined in Figure [5] 
r s and Tip are the meridional and circumferential stress, 
respectively, k v = sin ip/r is the circumferential curva¬ 
ture, and p — Ap g z is the hydrostatic pressure difference 
exerted on the capsule membrane. In this formulation, 
the bending stiffness has been neglected since the typical 
capsules used for these experiments are very thin and 
bendable [58] . 

In order to solve this system of ordinary differential 
equations numerically with a shooting method, the quan¬ 
tities A s and occurring on the right-hand side of the 
system must be calculated for given A v and r s by virtue 
of an elastic constitutive law. For the simple Hookean 
elastic law in eq [3] this was demonstrated in ref [58} 


A s 


- K ( T s - 7 ) - "(K - 1 ) + 1 , 

£/2D 

1^2 y s ((-V - 1 ) + v {K - 1 )) + 7- 


(14) 


The boundary conditions for the system of shape equa¬ 
tions are r(0) = 2 ( 0 ) = ip(0) = 0 and r(L 0 ) = a/2, where 
Lq is the contour length in the undeformed configuration 
and the end-point of the integration. The starting value 
t s ( 0 ) is free and serves as a shooting parameter to satisfy 
the boundary condition at L 0 [58 |. 

Here, we want to use the elasticity model for particle 
rafts developed above (eqs [3} [7] and [ 8 ]) , including the hard¬ 
core interactions. Therefore, depending on the size of 
compressive stresses, we have to switch between different 
constitutive elastic relations corresponding to Hookean 
constitutive relations or hard-core constitutive relations 
along the arc-length coordinate so as explained in detail 
in Appendix [B] 


A. Wrinkling and Modified Shape Equations 

In our model, we neglect the bending stiffness of the cap¬ 
sule membrane. The fact that the membrane is infinitely 
easy to bend implies that it will immediately wrinkle 
under compressive stresses. In ref 1581 it was shown that 
deflated pendant capsules exhibit wrinkles due to a com¬ 
pressive hoop stress t v < 0, as also indicated in Figure [5jx 
The capsule membrane is not axisymmetric in the wrin¬ 
kled region, but can be approximated by an axisymmetric 
pseudosurface around which the real membrane oscillates. 
The shape of this pseudosurface is determined by set¬ 
ting Tip = 0 where the original shape equations predict 
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negative hoop stresses. The assumption = 0 on the 
pseudosurface is common to various theories of wrinkling, 
for example tension field theory [ID] or far-from-threshold 
theory [|T]. This leads to a modified system of shape 
equations 

f'(so) = A s cos ip 
z'(s 0 ) = X s sin^> 

ip'{s 0 ) = ^-{p- Apgz) ( 15 ) 

T s 

r' s (s 0 ) = cos ip 

r 

where r is the radial coordinate of the pseudo-surface 
and t s the meridional stress measured per unit length of 
the pseudosurface, that is t s = TgX^/X^, when X v = f/r o 
denotes the pseudo hoop stretch. The system is closed by 
the equation 

A TsAy + E 2D - 7(1 + v) 

s E 2D - 2z/ 7 - 7 2 (1 - u 2 )/E 2 t> ‘ 1 J 

A discussion of these modified shape equations for the 
pseudo-surface and the automatic switching between nor¬ 
mal and modified shape equations during the integration 
can be found in ref f28l 

Wrinkling can also occur when the lattice is jammed on 
boundary 2 or point B. As in ref 1281 we handle wrinkling 
by introducing a pseudosurface (indicated with an overbar) 
and setting the hoop stress to zero. On boundary 2, we 
have X v = y^4L 2 /3 — X 2 /3 (note that X v refers to the 
hoop stretch of the real, wrinkled surface and is to be 
distinguished from the pseudo hoop stretch X v = f/ro). 
The condition t v = 0 is equivalent to 


O 

(A 

(A 


FIG. 6. Trajectories in the stress plane for integrations 
at different stages of deflation (see pictograms on the right). 
Thin dashed parts of the trajectories are jumps in r v . The 
integration always starts at the apex with r s — t v , that is, 
on the angle bisector, and runs through different domains of 
the stress plane. The parameters of the elastic model are 
A 2 D /7 = 5, v = 1/3, and L = 0.95. In the pictograms, thin 
horizontal lines indicate the transitions between the regions. 




Ts/l 


equation for r s . Note that the geometry of the pseudosur¬ 
face is not fixed by the condition X s = X v = L because 
the circumferential stretch X v of the pseudosurface is free. 
This is in contrast to the case of a lattice being stuck in 
point B without wrinkling as discussed above. 


B. Numerical Integration and Switching between 
the Shape Equations 


ri c) = -7 ’ 


E- 


2D 


1 — V 2 A< 


[(A v -l)+«/(A a -l)]. (17) 


We further have = r^ ■ X 2 /3X 2 . So the complete 

meridional tension, measured per unit length of the pseu¬ 
dosurface, is determined by 


f. = ^ 
A, 


7 (s) (A s ,A^ + ^r( c 


(18) 


This is a quite complicated function of A s , because r ^ 
and X v herein also depend on A s . It must be solved for X s 
to evaluate the right-hand side of the shape equations (eq 
151, which is done numerically in each integration step of 
the shape equations. 

When wrinkling occurs in point B, the shape equations 
for the pseudosurface (eq [l5]) can also be used, and since 
A s = X^ = L , the inversion of a stress-strain relation 
to obtain A s is not necessary. In point B, the hard-core 
stresses r| c ) and r ^ are independent of each other, and 
the spring contributions t| s ) and rI s ) are fixed because 
of X s = X v = L. From the wrinkling condition t v = 0 
we can, thus, calculate The meridional hard-core 
contribution r( c ) must be calculated from the differential 


The modified shape equations are integrated from the 
apex So = 0 to the attachment point So = L 0 to the 
capillary. On this way, the integration will run through 
different domains and must switch to the appropriate 
shape equations discussed in the previous section and 
Appendix [B] Figure [6] shows typical trajectories of the 
integration in the stress plane, that is, parametric plots 
of (r s (so),r ¥ ,(s 0 )) with s 0 £ [0,L o ]. 

We name the different domains of the stress plane as 
follows: The admissible domain (light green in Figure [6]) 
is abbreviated “A”, the red and blue ruled regions are 
termed “1” and “2” (because they stem from boundary 
1 and 2 in the stretch plane), and the white region is 
termed “B”. Regions 2 and B also appear in wrinkled 
form, and we then call them “2W” and “BW”. 

In the numerical integration, an event handler must be 
introduced which detects when the integration runs from 
one region into another. Changes from region A to regions 
1, 2 or B can be detected on the basis of strains, which 
are limited by the boundary of the forbidden domain as 
given by eq[6] The other direction, a change from a hard¬ 
core region into the A domain, occurs when the hard-core 
stresses become positive. 

Switching between the different hard-core regions is, 















however, a bit more complicated because the continuity 
conditions are less obvious. An elaborate variational 
calculation [42j shows surprisingly that t v may jump at 
transitions from B —> BW, B —> 2, 2 —>■ 1 and BW —> 1; see 
also Figure [Gj The physical reason behind this behavior 
is that t v is constitutively undetermined in region B: It 
may jump without requiring the hoop stretch X v to jump, 
which would be unphysical and would lead to a ruptured 
shape. This jump is necessary when starting in region B, 
where we see in Appendix[B]that the shooting parameter is 
eliminated: The jump, or rather its arc-length coordinate 
sj, serves as a substitute shooting parameter. In the 
following discussion of each trajectory shown in Figure [bj 
this becomes more clear. 

B — ► BW — > 1, Light Blue The transition from B to 
BW occurs at sj which can be chosen arbitrary 
(it just has to occur before the trajectory reaches 
the wrinkling region t v = 0). This is the shooting 
parameter that is adjusted to match the boundary 
condition t(Lq) = a/2 at the end of integration. 
The rest of the trajectory follows from the rules 
formulated above: The jump from BW to 1 occurs 
when the wrinkling condition X v < A v becomes 
false. As the stretches are fixed to L in B, BW and 
on the boundary of 1, the stretches are continuous 
at this transition, only the hoop stress jumps. 

B —> BW —> 2W —y 2 —y 1, Light Red.: The first tran¬ 
sition occurs again at a free position sj, and the 
remainder of the course follows: BW —> 2W is con¬ 
tinuous and occurs when the ratio r ^ /rj c ) becomes 
larger than 3; 2W —2 is also continuous and hap¬ 
pens when X v becomes larger than A^; and 2 —1, 
which has a jump in t v , is determined by X v be¬ 
coming larger than L. 

B —y 2 —y 1, Violet.: Again, the transition out of region 
B happens at a shooting parameter sj, and the sec¬ 
ond transition follows as explained in the previous 
trajectory. 

A — y 2 — y A, Dark Yellow.: For this trajectory, the 
shooting parameter is r s (0) as usual, because we are 
starting in region A. The integration switches from 
A to 2 when the boundary in the stretch plane is 
reached according to eq [6} and back to A when the 
hard-core stresses become positive. Both transitions 
are continuous. A transition from A to 1 is also 
possible and obeys the same reasoning. 

The dark green trajectory is trivial because it stays in re¬ 
gion A for the whole integration. More paths are possible 
and have been worked out, but the presented ones were 
most commonly met in the numerical investigations. 



FIG. 7. Wrinkle length as a function of the reduced volume 
for deflated shells with and without hard-core interactions, for 
v = 1/3 and E- 2 d as indicated in the legend. The shapes on 
the right for V = 0.7Vo reveal differences not only in L w but 
also in the overall shape (colors as in the legend, thin gray 
line is the undeformed shape, integration regions are indicated 
in the green contour). The five data points on the green line 
in the diagram indicate the shapes and trajectories shown in 
Figure [6] 


V. ANALYSIS OF COMPUTED DEFLATED 
SHAPES 

With the newly developed shape equations, we can 
compute deflated shapes of capsules according to the 
hard-core/soft-shell elasticity model. For the numerical 
analysis, we use 7 as the tension unit and a as the length 
unit. Dimensionless quantities are denoted with a tilde. 
Specifically, the shape equations depend on the reduced 
density difference p = a 2 Apg/'y , the reduced pressure 
p = ap/ 7 which controls the capsule volume and the 
reduced surface Young modulus i? 2 D = E^o/l- Their di¬ 
mensionless solutions contain the shape (r, z) = ( r/a , z/a) 
and tensions fj = Ti/j with i £ {s, p}. 

Starting with a Laplace-Young shape [25] with p = 0.25 
and po = 2, the pressure is lowered from p = po to p = 
O.lpo- Figure [7] shows the arc length L w of the wrinkled 
region as a function of the reduced volume for three such 
series of deflated shapes with different elastic parameters. 
The curves without hard-core interactions (red and blue, 
with L = 0) show that the onset of wrinkling occurs 
earlier for higher elastic moduli. If hard-core interactions 
are included and occur already for small compressions 
(green curve, L = 0.95), the wrinkling sets in early, even 
though the elastic modulus is small. Here, the wrinkling 
is induced by hard cores coming into contact. After hard¬ 
core contact, compressive stresses increase much more 
quickly with decreasing volume, as it can also be seen 
in the larger “stress loops” in regions 1, 2, and B in 
the trajectories shown in Figure [6] If these stress loops 
touch the gray shaded wrinkling regions, wrinkling is 
triggered (red and blue shapes in Figure [6]) . The quickly 
increasing stress loops after hard-core contact also give 
rise to the rather steep increase of the arc length L w of 
the wrinkled region as a function of the reduced volume 
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in the corresponding green line in Figure [TJ 

Thus, if one analyzes the shape of a deflated capsule, 
not knowing that it obeys the hard-core/soft-shell elas¬ 
ticity model but assuming that it is a usual Hookean 
membrane, the Hookean elastic modulus will be overesti¬ 
mated considerably. Figure [7] illustrates this, as the green 
line with A 2 d = 5 and L = 0.95 is much closer to the red 
line with A 2 d = 10 than to the blue line with E^d = 5. 


A. Shape Analysis of Theoretically Generated 
Shapes 

In refJ25] a shape analysis for deflated pendant capsules 
was developed. From experimental images, the contours 
are extracted and fitted with the solutions of the shape 
equations (eqs [13]). This allows one to determine the elas¬ 
tic modulus of the capsule membrane. The application of 
this method to bubbles coated with a layer of the protein 
hydrophobin revealed a very nonlinear elastic behavior, 
and the fitted elastic modulus as a function of the volume 
jumps at the onset of wrinkling. The essential features 
of the Hookean elastic modulus obtained in ref [25] are 
shown in Figure |8 Jd, where the area compression modulus 
A 2 d = A 2D /2(1 — v) is used instead of the surface Young 
modulus. 

We test whether this characteristic signature of the 
hydrophobin layer elasticity can be explained by our hard¬ 
core/soft-shell elasticity model. To this end, series of 
deflated shapes are computed using our modified shape 
equations for a given area compression modulus for 

the soft springs and a given hard-core length L. Then 
each theoretically generated shape is converted into a set 
of sampling points and fed to the usual shape analysis 
algorithm developed in ref [25] which uses a Hookean elas¬ 
ticity model without hard cores. The output of the shape 
analysis algorithm is a fitted Hookean area compression 
modulus A 2 d for each capsule volume, which is shown in 
Figure [8jil-a5 and demonstrates that the fitted Hookean 
area compression modulus A' 2 d can differ substantially 
from the value A or j g of the soft springs that has been used 
for the theoretical shape generation and that the fitted 
Hookean area compression modulus can change with vol¬ 
ume although A' or ig and the hard-core length L are fixed. 
In principle, the shape analysis algorithm using Hookean 
elasticity can also also give a fit value for the Poisson 
ratio v. Alternatively, we can fix v in the fit procedure 
to reduce the number of free parameters. 

The analysis is concentrated to soft shell area com¬ 
pression moduli of Aorig = 5 — 20 and hard-core lengths 
L = 0.95 — 0.99. Only for hard-core lengths L > 0.95, we 
will reproduce the jump in the fitted Hookean area com¬ 
pression modulus A" 2 d in the volume range V/Vq > 0.95, 
where it is also found for bubbles coated with a layer of 
the protein hydrophobin in ref l28l This motivates our 
choice of the parameter range for L. The range of hard¬ 
core lengths L > 0.95 is higher than the above estimates 
L ss 0.8 — 0.9 from the results of refs [Ml and 1551 for planar 


monolayers. These differences can be caused by different 
surface densities of hydrophobin proteins in the bubble 
geometry as compared to the planar geometry: In the 
preparation of bubbles [25] [37], hydrophobins can adsorb 
from the surrounding bulk solution to the interface, which 
might lead to an increased surface density (and, thus, an 
increased L) due to the additional adsorption energy as 
compared to the planar geometry [M] [50], where a fixed, 
amount of hydrophobin is spread onto the surface. In 
the fits using the Hookean elasticity model, we fix the 
Poisson ratio to the value v = 1/3, which is the appropri¬ 
ate value for a Hookean spring network, that is, for small 
compression if hard cores do not come into contact. 

When small hard-core lengths are combined with large 
area compression moduli, for example, L = 0.95 with 
A/mg = 20, we find that the hard-core interactions have 
no influence on the capsule shape, because the capsule 
starts to wrinkle before the hard-cores come into contact. 
The numerical integration then takes a path A —> AW —> A, 
and produces a shape that can also be produced by purely 
Hookean shape equations without hard-core interactions. 
Consequently, such shapes can be perfectly fitted with 
Hookean elasticity and the fitted Hookean area compres¬ 
sion modulus A2D agrees with the correct soft shell area 
compression modulus A or i g as expected (see, for example, 
L = 0.97 with A or ig = 15 in Figure [8^,3). In particular, 
there is no jump in the fitted Hookean modulus A' 2 d dur¬ 
ing deflation. Therefore, these cases are mostly omitted 
from Figure [8]rl-a5. There are some deflation series that 
are very close to the limit where the hard-core interactions 
cease to influence the shape (see L = 0.96 and A or i g = 10 
in Figure [8^2, for example), where there are only small 
deviations between fitted Hookean and original soft shell 
modulus. 

For computed deflation series in Figure [5[il a5 where 
the hard-core interactions profoundly influence the shape, 
the fitted Hookean area compression modulus indeed 
shows similar features as observed for hydrophobin cap¬ 
sules in ref [281 and shown in Figure §>• For small de¬ 
formations, where the hard cores are not yet in contact, 
the fitted Hookean modulus reproduces the original soft 
shell value, A 2D = A orig . When the hard cores come 
into contact, the fitted Hookean modulus K 2 d exhibits 
a jump and can grow much larger than the original soft 
shell value. For a better comparison between fit results 
for theoretical and experimental shapes, we show the di¬ 
mensionless ratio A 2 d / A or i g of fitted and soft shell area 
compression moduli in Figure]^) and c (for the experimen¬ 
tal fits in in Figure [8jo we identify A or i g with the mean 
of the fitted area compression modulus values before the 
jump). The peculiar dip in the fitted Hookean modulus 
A 2 d just before the jump is an artifact of fixing v = 1/3 
in the fits. Further tests with free v show that these 
points can also be fitted with a Hookean modulus A 2 d 
close to the original soft shell value A or i g but with the 
fitted v dropping to negative values; a result which lacks 
an intuitive explanation. Generally, the interpretation 
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FIG. 8. (al-a5) Fit results for theoretically generated capsule shapes of various area compression moduli 7\ or i g of the soft 

shells and hard-core lengths L. The fits are performed using the usual Hookean elasticity (without hard cores) and come to 
differing results K 211 for the fitted Hookean area compression modulus. In both shape generation and fits, the Poisson ratio is 
fixed to v = 1/3. Shapes exhibiting wrinkling are indicated by circles, unwrinkled shapes by squares, (b) Typical fit results for 
the dimensionless TGd / A or i g area compression modulus of a hydrophobin capsule as obtained in ref[28j The jump in the elastic 
modulus coincides with the onset of wrinkling, (c) Fit results to theoretically generated capsule shapes from subfigure (a4), 
normalized to 7\ or i g . (d) The relation between the plateau value A p i at of the fitted Hookean modulus TGd and the hard-core 
length L according to subfigures (al)-(a5). 


of the fitted value for v is not clear as it is obtained by 
using the inappropriate Hookean elasticity model also in 
the regime, where hard cores come into contact. The 
characteristic deflated volume, where the fitted Hookean 
modulus 7i2D exhibits a jump, increases with increas¬ 
ing hard-core length L but depends only weakly on the 
original soft shell modulus 7< or i g . For hard-core lengths 
L > 0.99 close to unity, unwrinkled shapes are practically 
unobservable. 

The deflation series in Figure [8}al — a5 also shows that 
contact of the hard cores triggers wrinkling (shapes ex¬ 
hibiting wrinkling are indicated by circles, unwrinkled 
shapes by squares). After the onset of wrinkling, the val¬ 
ues for the fitted Hookean modulus K 2 b reach a plateau 
and increase only slightly for increasing deflation in the 
computed shapes in Figure [8}\l-a5. For large L > 0.97, 
the plateau value is independent of the original soft shell 
modulus K OI i g . The plateau value strongly depends, 
however, on the hard-core length L as shown in Fig¬ 
ure §1- We can describe the observed plateau values 
by 7Fpi a t ~ 0.3(1 — L ) _1 . This can be explained by 
assuming that both the fitted Hookean elasticity and 
the hard-core/soft-shell elasticity have to describe the 


same sequence of shapes, the most prominent feature of 
which is the onset of wrinkles at a certain volume. If 
the onset of wrinkling is described by a fit with Hookean 
elasticity, the wrinkling criterion is t v = 0, which gives 
K 2 T>{X ip — 1) ~ — 7 , see eq[l4]for the hoop stretch at 
wrinkling. For hard-core/soft-shell elasticity wrinkling 
happens along the boundaries 1 or 2 according to the 
criterion from eq [6j which we approximate by X v ~ L. 
If Hookean elasticity is to describe the same shapes, 
X v as a function of the volume has to be identical, in 
particular, at the onset of wrinkling, which results in 
1 - Ap ~ 1/iFpiat ~ 1 - L or 7vpi a t ~ 1/(1 - L). Conse¬ 
quently, the fitted Hookean modulus K 2 d can exhibit a 
pronounced jump at the onset of wrinkling if the origi¬ 
nal soft shell modulus 7\ or i g is small and the hard-core 
length L is large; in Figure [8jr5 it jumps to more than 
its 6-fold value for A' orig = 5 and L = 0.99. The hard 
cores have an influence on the elastic properties if an 
increased plateau value is observed, that is, A' p i at > K or i g . 
Therefore, the domain of influence of hard cores is given 
by A'orig < 0.3(1 — A) -1 , which means sufficiently large 
hard-core lengths. 

We can use our findings for the plateau value of 
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the fitted Hookean modulues to extract the two hard¬ 
core/soft-shell elasticity parameters, the soft shell modu¬ 
lus A' or ; g = AT or ig /7 and the hard-core length L , from the 
Hookean fits of deflated shapes. First, the fitted Hookean 
area compression modulus at small deflation before wrin¬ 
kling can be used as an approximation to the soft shell 
modulus -Korig. For the hydrophbin capsule analyzed 
in ref B51 and shown in Figure [Sja, this gives a value of 
-Korig = 342mN/m. Second, we can use the above relation 
K^piat ~ 0.3(1 — A ) -1 to obtain the hard-core length A 
from the plateau values of the fitted Hookean modulus. 
The typical fit results for hydrophobin capsules as shown 
in Figure [ 8 jo do not show a plateau but a gradual decrease 
of the fitted modulus after the pronounced jump for large 
deflations and differ in this respect from the fits of hard¬ 
core/soft-shell capsules in Figure]^,. Extracting hard-core 
lengths A ss 1 — 0.3/A' p i at , from these fit results (using 
7 = 49.8mN/m), we find values of the hard-core length 
A which are decreasing from A = 0.99 at the jump to 
A = 0.98. This is a hint that the barrel in hydrophobin 
is not an ideal hard core but weakly compressible. 

Moreover, hydrophobin capsules also feature an ini¬ 
tial increase of the fitted Hookean modulus for small 
deflation (see Figure [8ju) which is absent in all the fits 
of hard-core/soft-shell capsules in Figure [8^1-a5. This 
could reflect an additional nonlinear stiffening of the soft 
hydrophobin shells upon compression. 

VI. CONCLUSIONS 

In this Article we developed an elasticity model for par¬ 
ticle rafts at an interface which consist of hard-core/soft- 
shell particles. Upon compression, the membranes formed 
by these rafts first behave according to the Hookean 
elastic law of the soft shells. Such “soft shell behavior” 
can be generated by any additional interaction between 
hard-core particles, in principle. If the hard cores come 
into contact, further compression is impeded. Additional 
stresses are then transmitted through the “skeleton” of 
hard cores, which must fulfill certain geometrical force 
balance constraints; see eq[ 8 j This strongly modifies the 
elastic response, that is, the constitutive stress-strain 
relations of the material, as soon as compression brings 
the hard cores of particles into contact. The model is 
characterized by an additional parameter A < 1, which is 
the ratio of hard core diameter to the equilibrium lattice 
constant. 

For a planar particle layer under compression, we find 
that this rather general elastic model of a particle raft 
gives rise to a stress-dependent apparent Poisson ratio, 
which increases upon compression beyond 1/3 and up to 
t'app = 1/(4A 2 — 1) if hard cores come into contact; see 
Figure [4] 

Curved interfacial layers are important for the elasticity 
of particle decorated capsules. We reviewed the shape 
equations for capsules attached to a capillary and modified 
them in order to include this hard-core/soft-shell elasticity 


model. This enabled us to compute deflated shapes. We 
find that the contact of hard cores by compression typi¬ 
cally triggers wrinkling of the capsule membrane because 
compressive stresses increase quickly with decreasing vol¬ 
ume after hard-core contact, as illustrated in the stress 
diagram in Figure [ 6 ] and the steep increase of the wrinkle 
length for decreasing volume in Figure [7] 

In a further analysis, theoretically generated shapes 
obeying the hard-core/soft-shell elasticity model were fit¬ 
ted with simple Hookean shape equations (without hard 
cores). The resulting fitted Hookean area compression 
modulus K 2 D, as shown in Figure (8}il-a5 for several pa¬ 
rameter combinations of soft shell modulus and hard-core 
length, exhibits features which are similar to the signa¬ 
tures obtained in ref[55]for hydrophobin capsules as shown 
in Figure [8]b. In particular we see a jump-like increase 
of the fitted Hookean modulus A' 2 d at a characteristic 
deflated volume and a plateau value of A 2 D for all smaller 
values. The jump in K 2 D coincides with the volume where 
hard cores start to come into contact, which also causes 
the onset of wrinkling. This explains the observed pro¬ 
nounced jump in the fitted Hookean area compression 
modulus at the onset of wrinkling with the effect of the 
hard cores. The characteristic deflated volume, where 
the jump in A^d occurs, increases with the hard-core 
length A approaching unity and is, thus, correlated with 
an increasing plateau value IUjd- This suggests that the 
Hookean model is only sufficient in a rather small volume 
range (above the jump in A' 2 d), where shapes are weakly 
deformed and unwrinkled. To describe wrinkled shapes 
adequately, we need to introduce the hard-core/soft-shcll 
elasticity characterized by the hard-core length A and 
the Hookean modulus AT or i g of the spring network. Nev¬ 
ertheless, the plateau value of AT 2 d characterizes the ef¬ 
fective Hookean compression resistance of the particle 
raft in situations, where hard cores come into contact. 
We demonstrated that values of the hard-core length A 
and the Hookean modulus A' or i g can already be inferred 
from a small set of distinct properties of the observed 
features of the fitted Hookean modulus A^d- For small 
deformations, before wrinkling and the associated jump 
of the modulus, we have AT or i g ss K^d- The value of A can 
be found from the plateau value of A 2 D after the jump. 

The hard-core/soft-shell model does not reproduce the 
initial increase of the fitted Hookean compression modulus, 
nor its final decrease as obtained in the hydrophobin fits. 
The initial increase of the fitted compression modulus 
might be reproducable if nonlinear spring interactions 
are included in the elasticity model (springs that stiffen 
upon compression). The decrease after the jump, on the 
other hand, might be considered as a decreasing hard¬ 
core length A, because for smaller A, the plateau value 
is smaller. This could be modeled by cores which can be 
slightly compressed. Replacing the hard-core interactions 
by spring-interactions with relatively large spring constant 
could achieve this. 

In summary, with respect to modeling the elastic prop¬ 
erties of hydrophobin coated interfaces, there are several 
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indicators that the proposed elastic model is too simple 
in distinguishing easily compressible domains and purely 
incompressible ones. A softer transition might produce 
results closer to the hydrophobin elasticity, modeled for 
example by springs whose spring constant slightly in¬ 
creases with compression, and then sharply increases to a 
large but finite hard-core spring constant. Such a model 
could even be theoretically more tractable because the 
resulting elastic law might involve a bijective mapping 
between stretches and stresses if the spring constant is a 
continuous function of the compression. 

The model is interesting not only with respect to hy¬ 
drophobin coated liquid interfaces but also as a generic 
model for rafts of interacting hard particles at liquid inter¬ 
faces, where the particle interactions give rise to the soft 
shell elasticity. Such type of membranes will exhibit a 
pronounced compression-stiffening after hard cores come 
into contact. For the capsules attached to capillaries 
this effect induced wrinkling and led to a corresponding 
jump in the apparent or fitted Hookean are compression 
modulus. 

The pronounced compression-stiffening of such mem¬ 
brane materials could also be used to stabilize structure, 
for example, closed capsules against compressive buck¬ 
ling [43] . The shape equations that we derived can be 
applied to the buckling behavior of particle decorated 
liquid droplets in future work as well. We expect that the 
interacting particle raft will give rise to high resistance 
to buckling if it is engineered such that hard cores come 
into contact at the critical buckling pressure. 
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Appendix A: Force-Equilibrium of the Hard-Core 
stresses 

We consider a jammed state of the lattice, that is, when 
( \ Xl X y ) is on the boundary of the admissible domain. 
Here, boundary 1 as specified by eq[6]and Figure [2}o is 
considered, where 3A 2 /4 + A 2 /4 = L. 

On this boundary, the lattice (see Figure [9^,) is com¬ 
pressed predominantly in the a:-direction. Not all hard 
cores of neighboring beads are in contact with each other, 
only those along the links that are drawn continuous in 
Figure HD 5 - The dashed link in this figure is a spring 
interaction, not a hard-core interaction, and is therefore 
ignored in the following. Figure[9]:, shows that the force F x 
applied to a bead is split into components F t tangential 
to the links. Trigonometric relations give cos a = F x /2F t . 
Analogous considerations give sin a = F y /2F t for the 
splitting of F y . Thus, we have 

p 

—T = tana (Al) 

F y 


as the condition that the skeleton is in static equilibrium. 

With Figure [9]i, we can relate the geometric angle a 
to the lengths of the links. The hard-core links have, by 
definition, length L. The vertical (dashed) link has a 
rest length of 1, and is stretched (or compressed) by the 
deformation to a length \ y • 1. We thus obtain 

tan a = —^ = X v _ ( A2 ) 

V^ 2 -^ Fy ^ L ' 2 -K 


Finally, we have to relate the forces F x and F y to the 
stresses and rff). Stresses are forces per length, and 
the investigated cell of the lattice has a height l y = Ay 
and width l x = l y /tana = V4L 2 — A 2 ; see Figure 
With rff) = F x /l x and = F y /l y , we thus arrive at 


7"( c ) A 2 

jy _ _ y 

r i c) 4 L 2 - A 2 ' 


(A3) 


With the strain constraint from eq[6]on boundary 1, this 
is equivalent to 


7"( c ) 1 A 2 

y _ __y_ 

3 \ 2 X ' 


(A4) 


Thus, only the ratio between the hard-core stresses is 
prescribed by the geometry of the lattice. 

Analogous results can be obtained on boundary 2, with 
all indices x and y interchanged. The complete result is 
therefore 


-(c) 

y 

-(c) 


A 2 /3A 2 



for A y > L 
for A v < L 


(AS) 


In point B with A x = A y — L , which lies on both 
boundaries 1 and 2, the lattice is uniformly compressed 
and all neighbouring beads are in contact. Equation (A5) 


then states that the ratio of the hard-core stresses is either 
3 or 1/3. In fact, due to the close packing of spheres, any 
value in between can also be realized, so that 


1 < T ’ 


3 t( c ) 


< 3 for \ T = A„ = L. 


(A6) 


Equations (A5) and (A6| give in the main text. 


Appendix B: Hard-core/soft-shell elasticity and 
capsule shape equations 

We want to use the elasticity model for particle rafts 
developed above (eqs [3} [7] and [8| , including the hard-core 
interactions, in the shape equations (eqs [l3]) for capsules. 
Therefore, depending on the size of compressive stresses, 
we have to switch between different constitutive elastic 
relations corresponding to Hookean constitutive relations 
(eqs[l4| or hard-core constitutive relations along the arc- 
length coordinate sq- 
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FIG. 9. Calculation of the ratio of the hard-core stresses 


We identify the ^-direction of the planar model with 
the (meridional) s-direction of the axisymmetric shell, 
and the y-direction with the (circumferential) (^-direction. 
With the help of the plots of the stretch and stress planes 
in Figure [3J: and d, we construct a suitable algorithm to 
calculate A s and t v from given X v and t s : 

(i) We check if the point is in the admissible domain 
or on the boundary by calculating A^ b )(A v ) with 
eq Kh and (A^, A v ) using eq kl] which is the 
smallest possible stress in the admissible domain. If 
the given r s is larger (less compressive) than this 
value, the point is in the admissible domain, if it is 
smaller (more compressive) then the point is on the 
boundary to the forbidden domain. 


equations. The reason is that the confinement to A s = 
A v = L already determines the shape of the capsule: It is 
uniformly compressed. The circumferential stretch r/ro = 
L directly implies r(so) = Tro(so)- Inserting this solution 
into the differential equation for r in the system of shape 
equations (eqs 13) then yields i/j(sq) = ipo(so), where 


ipo is the slope angle in the undeformed configuration. 
The differential equation for z then becomes ^(so) = 
Lsin^o and has the solution z(so) = Lzo(sq) + c with 
some constant c depending on the starting value. From 
the differential equation for ij} in the shape equations (eqs 
131 we can then deduce 


= + p- Apg z'j /k v 


(Bl) 


(ii) If the point is in the admissible domain, eqs 14 


be used to calculate A s and t v . 


can 


(iii) If the point is on the boundary, then we know that 
A s = A( b i(A v ) according to eq[bj In addition, we 
can calculate the spring contributions (A^ b \ A v ) 
from Hooke’s law (eq[3]). The hard-core contribution 
can then be obtained as the difference between the 
given r s and the spring contribution, r| c ) = t s — 
r| s ). This value should be negative. The sought 
stress can then be calculated from the Hookean 
contribution according to eq [3l and the hard¬ 
core contribution r ^ according to me force balance 
condition (eq[8j), so that t v = (X^\X v ) + 

Thus, we can use the above shape equations (eqs 13]) also 
for the computation of deformed shapes for capsules obey¬ 
ing our hard-core/soft-shell elasticity model. During the 
integration of the shape equations, however, the closing 
relations (eqs[l4]) that are necessary to compute the right- 
hand side must be replaced by the above procedure (iii) 
in regions where the lattice is on boundary 1 or 2. 

This method does not work when the lattice is stuck in 
point B, which must be handled separately in the shape 


by inserting the known solutions, where k So is the merid¬ 
ional curvature in the undeformed configuration. So in 
principle, only t s must be determined by solving its differ¬ 
ential equation. In order to keep the code of the numerical 
implementation consistent, however, we solve the full sys¬ 
tem of shape equations (eqs [l3| with the closing relat ions 
(eqs[l4]) modified to contain the explicit result, eq Bl just 
derived. This produces the correct numerical solutions in 
the same “data format” as in all other parts. 

When the lattice is in the jammed state B already at 
the start point of integration (at the apex), we need to 
evaluate the limits of the explicit solutions for sq —> 0. At 
the apex, the meridional and circumferential curvatures 
coincide, k s (0) = /^(O) = k S o /L = p 0 /2~/L, where the 
last step can be derived from the Laplace-Young equation 
and po is the pressure inside the capsule in its undeformed 
configuration. From the force balance, the meridional and 
circumferential tensions follow as r s (0) = r v (0) = pLj/p 0 . 
This finding has a remarkable impact: t s (0) is fixed by 
the external parameters, and cannot serve as a shooting 
parameter. The problem of having lost the only shooting 
parameter is resolved in the main text in the discussion 
of continuity conditions. 
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